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ABSTRACT 

We present new Gemini-North/NIRI and VLT/SINFONI //-band spectroscopy for a flux limited sample of 
40 z ~ 4.8 active galactic nuclei, selected from the Sloan Digital Sky Survey. The sample probably contains 
the most massive active black holes (BHs) at this redshift and spans a broad range in bolometric luminosity, 
2.7 X 10"*^ < Lboi < 2.4 X lO'^^ergss"^ The high-quaUty observations and the accurate fitting of the Mg II A2798 
line, enable us to study, systematically, the distribution of BH mass (Mbh) and normalized accretion rate 
(L/LEdd) at z ~ 4.8. We find that 10** < Mbh ^ 6.6 x 10^ with a median of - 8.4 x 10** M©. We also find 
that 0.2 ;S L/LEdd ^ 3.9 with a median of ^ 0.6. Most of these sources had enough time to grow to their 
observed mass at z ~ 4.8 from z = 20, assuming a range of seed BH masses, with ^ 40% that are small enough 
to be stellar remnants. Compared to previously studied samples at z ~ 2.4 and ~ 3.3, the masses of the z — 4.8 
BHs are typically lower by ^ 0.5 dex. and their L/Ledd is higher by a similar factor. The new z — 4.8 sample can 
be considered as the progenitor population of the most massive BHs at z — 2.4 and ~ 3.3. Such an evolutionary 
interpretation requires that the growth of the BHs from z — 4.8 to z — 3.3 and z — 2.4 proceeds with short duty 
cycles, of about 10-20%, depending on the particular growth scenario. 
Subject headings: galaxies: active - galaxies: nuclei - quasars: emission lines 



1. INTRODUCTION 

The local Universe provides ample evidence for the exis- 
tence of Super-Massive Black Holes (SMBHs) in the centers 
of most galaxies. Typical masses are in the range Mbh ^ 
10^- 10^ Mq, with few exceptionally massive objects reach- 
ing lO'^Mfj. As first argued by Soltan (1982), the total lo- 
cal BH mass is consistent with the total radiation emitted by 
accreting SMBHs in active galactic nuclei (AGNs). The ac- 
cumulation of mass onto SMBHs can be traced back through 
cosmic history by analyzing the redshift-dependent quasar lu- 
minosity function. Such studies suggest that the peak epoch 
of SMBH growth was at z 2-3 (e.g., Miyaji et al. 2001; 
Hasinger et al. 2005; Silverman et al. 2008; Croom et al. 
2009). However, this statistical approach does not provide 
sufficient information about the mass of individual SMBHs 
at various redshifts. A more detailed evolutionary study re- 
quires such measurements, in combination with a reliable es- 
timate of the bolometric luminosity (Lboi) and hence the light- 
to-mass ratio or, equivalently, the normalized accretion rate 

L/ LEdd = Lhol/ ^Edd- 

Several recent studies suggest that the more massive BHs 
experience most of their growth at very early epochs (z ^ 3), 
while those observed as AGNs in the the local Universe tend 
to have lower Mbh ("downsizing"; e.g., Marconi et al. 2004; 
Shankar et al. 2009). A similar effect is suggested for the dis- 
tributions of L/Ledd, such that at z^i2 most AGNs accrete 
close to their Eddington limit (Merloni 2004; Shankar et al. 
2009). It is also predicted that there should be an anticorre- 
lation between Mbh and L/LEdd at all redshifts. These trends 
are in general agreement with observations (e.g., McLure & 
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> 10 Mf^,), to account for 
6 (e.g.. Fan et al. 2006; 



Dunlop 2004; Netzer & Trakhtenbrot 2007; Shen et al. 2008). 
Detailed simulations of galaxy mergers predict that the in- 
stantaneous AGN luminosity and SMBH accretion rate vary 
greatly on very short timescales (Di Matteo et al. 2005; Hop- 
kins et al. 2006; Sijacki et al. 2007). The overall BH accretion 
period following major mergers lasts ^ 1 Gyr, out of which 
the central source would appear as a luminous, unobscured 
AGN for at most a few x lOOMyr. Mbh may grow by as 
much as a factor ^ 1,000 during such mergers. The merger 
history of SMBH hosts can also be traced in cosmological 
simulations of structure formation, by following halo merger 
trees (e.g., Volonteri et 2003). Almost all these models re- 
quire (or assume) that at z 3 most AGN would accrete close 
to, or indeed at their Eddington limit. The fast growth has to 
last almost continuously from very early epochs (z ^ 20) and 
involve massive seed BHs (Mgeed 
the very massive BHs detected at z 
Volonteri 2010, and references therein). The combination of 
structure formation models predictions at z ^ 3 - 6 with the 
recently observed high clustering of high-redshift luminous 
AGNs (Shen et al. 2007) suggests that the typical duty cy- 
cles - the fraction of the total time involving fast accretion - 
should remain above 0.5 and probably reach unity (e.g. White 
et al. 2008; Wyithe & Loeb 2009; Shen et al. 2010; Shankar et 
al. 2010a; Bonoli et al. 2010). These studies use various pre- 
scriptions to link AGNs with their dark matter halos. At lower 
redshifts (i.e. z ^ 2), the accretion can be more episodic with 
a typical duty cycle between ^ 10~^ and ^0.1 (e.g., Marconi 
et al. 2004; Merloni 2004; Shankar et al. 2009) A compre- 
hensive, up-to-date review of many of these issues is given in 
Shankar (2009). 

In order to test these scenarios, Mbh and L/Lgdd ought 
to be measured in large, representative samples. For unob- 
scured, type-I AGNs, this is usually achieved by using "single 
epoch" (or "virial") Mbh determination methods, which are 
based on the results of long-term reverberation mapping cam- 
paigns. These methods are based on estimating the size of the 



2 



TRAKHTENBROT ET AL. 



broad line region (BLR) and involve an empirical relation of 
the form Rblr °= (ALa)", where XL\ is the monochromatic 
luminosity in a certain waveband. Combining Rblr with 
the assumption of virialized motion of the BLR gas, we get 
Mbh = /G^~^^"^BLR' where / is a geometrical factor of order 
unity (e.g., Kaspi et al. 2000; Vestergaard & Peterson 2006; 
Bentz et al. 2009). The effect of radiation pressure force on 
such estimates is still a matter of some discussion (Marconi 
et al. 2008), but recent work suggests it is not very impor- 
tant (Netzer 2009a; Netzer & Marziani 2010). Single epoch 
mass estimate methods based on the H/? and Mg II A2798 
lines (e.g., Kaspi et al. 2005 and McLure & Dunlop 2004, 
respectively) were used to estimate Mbh up to z 2± 2 in large 
optical surveys (e.g., Corbett et al. 2003; McLure & Dunlop 
2004; Netzer & Trakhtenbrot 2007; Fine et al. 2008; Shen 
et al. 2008). Much smaller samples of z > 2 sources were 
studied by observing H/3 or Mg II in one of the NIR bands 
(Shemmer et al. 2004, hereafter S04; Kurk et al. 2007, here- 
after K07; Netzer et al. 2007, hereafter N07; Marziani et al. 
2009; WiUott et al. 2010, hereafter WIO). 

Mbh can also be estimated from the broad C IV A 1549 line, 
using specifically calibrated relations (e.g., Vestergaard & Pe- 
terson 2006). This would potentially enable the study of 
large samples of AGN at high redshifts. However, there is 
clear evidence that C iv-based estimates of Mbh are unreli- 
able. Baskin & Laor (2005) found that the C IV line is often 
blue-shifted with respect to the AGN rest-frame, which sug- 
gests that the dynamics of the C iv-emitting gas may be dom- 
inated by non-virial motion. Several studies of large samples 
clearly demonstrate that the relation between the widths of the 
C IV line and of the lower ionization lines (H/3 and Mg ll) is 
weak and shows considerable scatter (e.g., Shen et al. 2008; 
Fine et al. 2010) that is inconsistent with the virial assumption 
used in such mass estimators. Finally, our own study (N07) 
of luminous z — 2.4 and ~ 3.3 AGNs shows a complete lack 
of correlation between C iv-based and H/3-based estimates of 
Mbh- These discrepancies become crucial at high redshifts 
and large Mbh and lead to the conclusion that only H/3-based 
and Mg ll-based mass estimates are reliable enough to infer 
the properties of such sources. 

Our previous project (S04 & N07) presented the largest 
sample of z ~ 2.4 and ~ 3.3 type-I AGNs for which Mbh and 
L/Ledd were reliably measured using NIR H/3 spectroscopy. 
The distribution of L/Ledd at z ~ 2.4 and ~ 3.3 was found 
to be broad, and about half of the sources had L/L^fi^ < 0.2, 
inconsistent with several of the models mentioned above. In 
particular, the typically low accretion rates and the very high 
masses (up to log (Mqh/Mq) ~ 10.5) also mean that ^ 60% 
of the sources did not have enough time to grow to the ob- 
served Mbh by continuous accretion at the observed rates. 
These findings suggest that an epoch of faster SMBH growth 
must have occurred, for most objects, at z > 3.5. In order to 
probe such redshifts, the Mg II line must be observed in either 
the H or the /T-bands. Practically, this corresponds to focus- 
ing on z — 4.8 or z — 6.2 sources. 

In this paper we present a systematic study of Mbh and 
L/Ledd in a large, well-defined sample of z — 4.8 type-I 
AGNs. This is based on new //-band spectroscopic obser- 
vations, which enable the measurement of the Mg II line. We 
describe the sample selection and the observations in ^and 
the way we deduced Mbh and L/LEdd in ^ The main results 
of these measurements are presented in ^and discussed in ^ 
where we compare these results to those of other high-redshift 



samples. The main findings are summarized in ^ Through- 
out this work we assume a standard ACDM cosmology with 
l^A = 0.7, = 0.3 and //o = 70 km s'^ Mpc'^ . 

2. SAMPLE SELECTION, OBSERVATIONS AND DATA REDUCTION 
2.1. Sample selection 

Sources were selected to allow measurement of the 
Mg II A2798 emission line and the continuum flux at 3000A 
(^300o), thus providing reliable Mbh and L/Lndd estimates. 
The sources were selected from the sixth data release (DR6; 
Adelman-McCarthy et al. 2008) of the Sloan Digital Sky Sur- 
vey (SDSS; York et al. 2000). We first limited our search to 
"QSO"-class objects with z 4.65-4.95, so that the Mgll 
line and /3000 could be observed within the //-band. The 
completeness rate of SDSS spectroscopy for QSOs at this red- 
shift is above 98%, down to a limiting magnitude of / 2± 20 
(Richards et al. 2006a). This initial search resulted in 177 ob- 
jects. Next, we verified that no broad absorption features are 
present near the SDSS-observed CiV lines, since these are 
also expected to appear in the Mg II profile. This reduced the 
list of candidates to 129 sources. 

To estimate the //-band flux of the targets, we ex- 
trapolated the rest-frame flux density at 1450A to the 
observed-fiume flux density at 1.65/im, using the fiy 
i/""^^ spectral energy distribution (SED) presented by Van- 
den Berk et al. (2001). We applied a flux limit of 
/a (1450) ^6 X 10~'^ergscm~^s~^~', which translates to 
fx{l.65fim) ^2.1 X 10~'^ergscm~^s~' ~^ to include only 
sources which would provide NIR spectra with 5'/A^^5-10 
within reasonable exposure times. This formal flux limit 
translates to Lboi ~ 5 x 10^^ ergs s~' (see below). Our flux 
limit results in the omission of only 6 candidates, leaving 123 
sources. Different observational constraints (a combination 
of declination, observational seasons etc.; see § 32.2l i. forced 
us to observe only 40 of these 123 candidates. We verified 
that the distribution of f\ (1450) for the observed sample is 
very similar to that of the flux-limited sample of 123 sources, 
as well as the entire initial sample of 177 SDSS objects. In 
particular, we find 4 sources (^ 1 1 % of the observed sample) 
with /a (1450) ^7.5 x 10~'^ergscm~^s~'~' and 15 sources 
(-43%) with /a (1450) < 1.5 x lO-^ergscm-^s"' "i. Given 
this and the high level of completeness in the SDSS, we con- 
sider this sample to be representative of the complete sample 
of z ~ 4.8 luminous AGNs. 

Several studies suggest that the mass and accretion rate evo- 
lution of SMBHs may be connected with their radio proper- 
ties (see, e.g., McLure & Jarvis 2004; Shankar et al. 2010b, 
but also Woo & Urry 2002). To determine the radio properties 
of the z — 4.8 sources, we utilized the cross-matched catalog 
Faint Images of the Radio Sky at Twenty-Centimeters radio 
survey (FIRST; Becker et al. 1995; White et al. 1997). Due 
to the high redshift of our sample, and the relatively low sen- 
sitivity of FIRST, most candidates have only upper limits on 
their radio fluxes, and thus upper limits on the radio-loudness 

(/?L = YJ^^'^ Kellermann et al. 1989). For those sources, 
we estimated /?l based on a 3(7 FIRST upper flux limit of 
fij (5GHz) = 0.6 mJy and assuming a radio SED of 0= i/"" **. 
fi, (4400) was estimated from fx (1450) and the Vanden Berk 
et (2001) template. All but five observed targets (see Ta- 
ble[T]) have upper limits on /?l which are below — 40. Two of 
the remaining sources (J0210-0018 & J1235-0003) have firm 
FIRST detections (f^ [5GHz] = 9.75 and 18.35 mJy, respec- 
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Table 1 

Observation Log 



Object ID (SDSS J) 


^SDSS 


111 s tru mc nt 


Obs. Date 


Tot3.1 Exp 


//-band magnitude 










Time (sec) 


b 

spectro. 


imaging' 


000749.17+004119.4 


4.837 


SINFONI 


6 & 29 Jun. 2008 


4200 


18.59 


n8 1^ 


003525.28+004002.8 


4.757 


SINFONI 


29 & 30 lun. 2008 


4200 


17.82 


(18 0) 


021043.15-001818.2 


4.733 


NIRI 


4 Sep. 2007 


6380 


18.73 


n7 8"! 


033119.67-074143.1 


4.738 


SINFONI 


1 Oct. 2008 


2400 


17.49 




075907.58+180054.7 


4.861 


NIRI 


18 Oct. & 11 Nov. 2007 


8100 


17.31 




080023.03+305100.0 


4.687 


NIRI 


7 Dec. 2008 


6160 


17.06 




080715.12+132804.8 


4.874 


NIRI 


21 Dec. 2008 


4590 


17.52 




083920.53+352457.6 


4.777 


NIRI 


15 Dec. 2008 


7290 


18.26 




085707.94+321032.0 


4.776 


NIRI 


6 Ian. 2008 


6490 


17.06 




092303.53+024739.5 


4.660 


SINFONI 


23 Dec. 2008 & 1 Jan. 2009 


6000 


18.44 


(18 4) 


093508.50+080114.5 


4.699 


SINFONI 


20 Dec. 2008 


2400 


17.89 


17.9 


093523.32+411518.7 


4.836 


NIRI 


7 & 12 Dec. 2008 


6440 


17.36 




094409.52+100656.7 


4.748 


SINFONI 


23 Mar. 2010 


2400 


17.77 




101759.64+032740.0 


4.917 


SINFONI 


28 Nov. 2008 & 2 Jan. 2009 


6000 


18.63 


fl8 91 


105919.22+023428.8 


4.735 


SINFONI 


17 Feb. 2009 


2400 


17.75 


18.0 


111358.32+025333.6 


4.882 


SINFONI 


6 Apr. 2008 


3000 


17.9 




1 14448.54+055709.8 


4.793 


SINFONI 


6, 8 & 22 Apr. 2008 


8100 


18.56 


18.6 


115158.25+030341.7 


4.698 


SINFONI 


22 Apr. 2008 


6000 


18.94 




120256.44+072038.9 


4.785 


SINFONI 


3 & 6 Ian. 2009 


4800 


18.15 


18.0 


123503.04-000331.6 


4.723 


SINFONI 


4 May. 2008 


8700 


18.46 




130619.38+023658.9 


4.852 


SINFONI 


6 Feb. 2009 


3300 


16.81 


16.9 


131737.28+110533.1 


4.810 


SINFONI 


20&21 Ian. 2009 


4500 


18.01 


ri7 91 


132110.82+003821.7 


4.716 


SINFONI 


7 & 9 Apr. 2008 


3000 


18.4 




132853.67-022441.7 


4.695 


SINFONI 


8 Feb. 2009 


2400 


18.07 




133125.57+025535.6 


4.737 


SINFONI 


14 & 26 Apr. 2008 


4800 


18.67 


18.7 


134134.20+014157.8 


4.670 


SINFONI 


8 Feb. 2009 


1800 


17.02 




134546.97-015940.3 


4.714 


SINFONI 


26 Apr. 2008 


1800 


18.38 




140404.64+031404.0 


4.870 


SINFONI 


28 Mar. 2010 


4500 


17.71 


17.7 


143352.21+022714.1 


4.721 


SINFONI 


29 Mar. 2010 


600 


16.7 




143629.94+063508.0 


4.850 


SINFONI 


28 Mar. 2010 


2400 


17.79 


18.0 


144352.95+060533.1 


4.879 


SINFONI 


4 Apr. 2010 


4800 


18.56 


18.9 


144734.10+102513.2 


4.686 


SINFONI 


7 & 8 Apr. 2010 


5400 


18.78 


(18.6) 


151155.98+040803.0 


4.686 


SINFONI 


4 Apr. 2010 


2400 


17.99 


17.8 


161622.11+050127.7 


4.872 


SINFONI 


8 Apr. 2010 


900 


16.89 


17.3 


165436.86+222733.7 


4.678 


NIRI 


19 Aug. 2007 


3480 


17.66 




205724.15-003018.0 


4.663 


NIRI 


12 Aug. 2008 


5800 


16.77 




220008.66+001744.8 


4.818 


NIRI 


12 Oct. 2007 


6160 


17.51 


(17.5) 


221705.72-001307.7 


4.689 


SINFONI 


19 Apr. 2008 


3300 


18.12 


(18.4) 


222509.16-001406.8 


4.888 


NIRI 


3 Sep. 2007 


6380 


17.15 


(17.1) 


224453.06+134631.8 


4.657 


SINFONI 


6 & 8 lul. 2010 


4800 


18.67 





" Redshift obtained from the SDSS archive, based on rest-frame UV emission lines, 
//-band (Vega) magnitude derived from the calibrated spectra by synthetic photometry. 

//-band (Vega) magnitude derived from direct CTIO/ISPI imaging, or from the UKIDSS database (in parentheses). 



lively) and are thus considered as radio loud AGNs — 104 
and 1080, respectively). Three additional targets were not ob- 
served by FIRST and thus have no viable radio data. In what 
follows, we verified that the inclusion or exclusion of these 
5 sources does not significantly affect our results. In sum- 
mary, the 40 z ~ 4.8 AGNs presented here comprise a flux 
limited sample which represents a large fraction of all such 
SDSS sources. In terms of bolometric luminosity, it is com- 
plete down to Lboi — 5 X lO'*^ ergs s~' and most of the sources 
are not radio-loud. The basic properties of the 40 sources are 
given in Table [T] 

2.2. Spectroscopic observations and reduction 

The //-band spectra of our z — 4.8 sample were obtained 
using the Gemini-North and the Very Large Telescope (VLT) 
observatories. The campaign was split over several semesters, 
with brighter targets predominantly observed with the Near 
Infra-Red Imager instrument (NIRI; Hodapp et al. 2003) on 
Gemini-North (as part of programs GN-2007B-Q-56 and GN- 
2008B-Q-75) and fainter ones with the more sensitive SIN- 
FONI instrument (Eisenhauer et al. 2003) on the VLT-UT4 (as 
part of programs 081.B-0549, 082.B-0520 and 085.B-0863). 



The log of observations is given in Table [T] All the obser- 
vations were performed in queue/service modes, requested 
not to exceed airmasses of ~ 1.7 and seeing of ^ 1", dur- 
ing clear skies and "gray-time" lunar phase. The 1 1 Gemini- 
N/NIRI targets were observed through a 0.75" x 1 10" slit and 
the G5203 grism at the f/6 setup, providing R ^ 520. The 
sub-integrations were of ^ 300 sec, with dithers of 6"along 
the slit after each sub-integration. The 29 VLT/SINFONI 
targets were observed through the 8" x 8" FOV (a.k.a "250 
mas/spaxel") We used SINFONI's "Hh-K" mode, since the 
broad Mg II A2798 line can be well-resolved with the result- 
ing R 1500. The /T-band spectra assisted in determining 
the continuum flux of the SINFONI targets. Sub-integrations 
were of 150 sec. or 300 sec. and (diagonally) dithered by 
^ 6" across the FOV. Despite this type of dithering, the spec- 
troscopically resolved OH sky emission contributes a major 
source of noise to our reduced data (see below). In both ob- 
servatories, telluric standards of spectroscopic types B, A and 
G, chosen to have similar airmasses, were observed immedi- 
ately before or after the science targets. 

The reduction of the raw data was carried out using the 
standard pipelines of the respective facilities. We used the 
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Figure 1. Rest-frame spectra of the z ~ 4.8 sources under study. The spectra are ordered by increasing right ascension. In each panel, we show the observed flux 
density (thin blue Hne) and the best-fit model (thick black line), which is composed of a continuum component (dotted line), an Fe II & Fe III emission complex 
(thin dashed curve), and the total (two-component) Mg II A2798 line (thin black lines). 



GEMINI V. 1.8 package in IRAF to co-align and combine 
the sub-integrated frames while subtracting the sky emission. 
The one-dimensional spectra were extracted through a typi- 
cal aperture of 20 pixels which correspond to ^ 2.3". The 
SINFONI V.2.0.5 pipHne (under the ESOREX environment) 
was used to extract "3D data-cubes" from each sub-integrated 
frame, which include an individual wavelength calibration. 
The data-cubes were then co-aligned and combined, while 
the sky emission was estimated from the appropriate dithered 
pointings. The combined data-cube of each source was ex- 
amined to estimate the actual PSF, and an appropriate circular 
aperture was used to extract the 1-d spectra. Due to the vary- 



ing flexture of the instrument during the long (Ihr) SINFONI 
Observing Blocks, there is a noticeable "wavelength flexture" 
effect in the extracted 1-d spectra. This effect results in P- 
Cygni-like spectral features and can be treated separately (see 
Davies 2007). We note that the observable properties critical 
to this work (i.e. L3000 and FWHM[Mg ll]) are not very sen- 
sitive to such narrow ( < 400kms-i)OH -originated features. 

The spectra of the standard stars were reduced and extracted 
in the same manner as above, using the same apertures as 
those used for the science targets. The published spectral 
types of these stars were used to correct for the instrumen- 
tal response of each science observations and their published 
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Figure 1. -continued. 



(2MASS; Cutri et al. 2003) magnitudes were used to flux- 
calibrate the science spectra. The typical photometric error 
associated with the 2MASS magnitudes is ^ 10%. 

2.3. Photometric observations 

Wide-field H and /T-band imaging for 11 of the z — 4.8 
sources was obtained at the CTIO Blanco telescope, with the 
ISPI instrument (van der Bliek et al. 2004). We aimed at hav- 
ing at least 3 bright 2MASS stars observed in the same field 
each of the science targets. Each field was imaged to achieve 
S/N ^ 20 for the science targets, through a series of dithered 
short exposures, to avoid saturation. The CTIO observations 
were carried out with mostly clear conditions, with seeing of 
~ 1", during a single run on the night of March 4th., 2010. 



The raw imaging data was reduced following standard proce- 
dures, including bad pixel removal, dark and flat calibration 
and co-alignment of the sub-exposures. Aperture photometry 
was performed over the science targets and the several visi- 
ble 2MASS stars in each field. Whenever possible, we used 
only stars with best-quality 2MASS photometry for the rela- 
tive scaling of the science target fluxes. The magnitudes thus 
obtained are also listed in Table[Tl for the relevant sources. 

Since the spectroscopy and imaging observations were not 
simultaneous, there exists a real uncertainty regarding pos- 
sible flux variability of our z — 4.8 sources. Our targets 
are relatively luminous, and thus are not expected to have 
large variability amplitudes (e.g. Bauer et al. 2009 and refer- 
ences therein). We compared the photometric and the spectro- 
scopic fluxes, achieved by synthetic photometry of the SIN- 
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Figure 1. -continued. 

FONI & NIRI spectra. The median discrepancy is ;S 0.09 
mags, in the sense that the (later observed) photometric fluxes 
are on average fainter than the spectroscopy-deduced fluxes. 
The standard deviation is close to 0.2 mags. Photometric 
fluxes for 10 additional sources were obtained through the 
fourth data release of the UKIRT Infrared Digital Sky Sur- 
vey (UKIDSS/DR4; Lawrence et al. 2007). The photome- 
try of these sources is also consistent with the fluxes deduced 
from our spectrophotometric calibration. Excluding one sig- 
nificant outlier (J02 10-00 18), the mean deviation between the 
two methods is, again, close to 0.2 mags. 

With this evidence of little flux variation, and reliable spec- 
trophotometry for a large fraction 50%) of arbitrary se- 
lected sources, we conclude that the accuracy of our spec- 
trophotometric flux calibration is high, with uncertainties of 
about 10-15%. In order to correctly probe the instantaneous 
emission from the sources, we thus adopt the spectrophoto- 
metric fluxes to measure luminosities, even for those sources 
which were photometrically observed. The full set of cali- 
brated spectra is presented in Figure[T] 

3. LINE FITTING AND Mbh & i/i-Edd DETERMINATION 

3.1. Line fitting 

To determine Mbh and L/L-^iii we fit the observed spec- 
tra with a model which combines a continuum component, 
Fe Hand Fe III lines, and the two doublet Mg II lines. The 
procedure used here is similar to the one presented in several 
other studies (e.g. Shen et al. 2008 and references therein). 
It is based on a new code that was developed and extensively 
tested on a lai-ge sample of 0.5 ;S z ^ 2 SDSS type-I AGNs, 
and will be described in detail in a future publication (Trakht- 
enbrot et al. 2011). In particular, this large sample was used to 
verify that important quantities such as narrow emission lines 
and bolometric corrections, are consistent with those derived 




2600 2700 2800 2900 3000 
Rest-frame wavelength [A] 



from the well-studied H/3-[0 III] A5007 emission region (rest- 
frame wavelengths of 4600-5100A). This was done using a 
sub-sample of ^ 5,000 sources at 0.5 ;S z ^ 0.75, where the 
SDSS spectra show both the and Mg II fines. 

The fitting of the z ~ 4.8 AGNs was preformed individu- 
ally, verifying a satisfactory match between the data and the 
model. As a first step, a linear pseudo-continuum is fit to the 
flux around 2655 and 3020A. Next, we fit the Fe II and Fe III 
emission complexes redward and blueward of Mg II. For this 
we use a template made of the composite prepared by Vester- 
gaard & Wilkes (2001) that was supplemented by several pre- 
dicted Fe II lines (kindly provided by G. Ferland). These are 
mainly Fe II emission lines which coincide in wavelength with 
the Mg II line, where the observationally-based template of 
Vestergaard & Wilkes (2001) is incomplete. This part of the 
template is very similar to the models presented in Sigut & 
Pradhan (2003; their Fig. 13) and in Baldwin et al. (2004; 
their Fig. 5). In both cases, the Fe II flux under the Mg II line 
is dominated by the red extension of the ~ 2750A emission 
complex. The additional Fe II flux tends to flatten in sources 
with very broad lines, but the overall effect on the Mg II line 
fitting is very small. We consider this template to be more re- 
liable than the addition of constant flux under the Mg II line, 
adopted by Kurk et al. (2007) and Fine et al. (2008). Sev- 
eral other studies, such as Salviander et al. (2007) and Shen 
et al. (2008), use templates very similar to ours, following the 
models in Sigut & Pradhan (2003). 

Our basic iron template is based on a single line profile 
with FWHM ~ 1 , 175 kms~' . To account for the range of ob- 
served line widths, we created a grid of broadended Fe tem- 
plates, by convolving the basic template with single Gaussian 
profiles of varying width, ranging from FWHM = 1,200 to 
lOjOOOkms"'. We note that the templates that correspond to 
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FWHM > 4, 000km s"^ lack almost any distinguishable emis- 
sion features. The grid of broadened Fe templates is fitted to 
the continuum-subtracted spectra over the rest-frame wave- 
length regions of 2600-2700A and 2900-3030A, and the 
best-fit template is chosen by standard minimization. We 
allow the Fe template to be shifted with respect to the systemic 
redshift and for its width to differ from that of the Mg II line. 
After the subtraction of the best-fit Fe template, the pseudo- 
continuum is re-fitted and we iterate the Fe-fitting process 
once again. This is done to ensure the convergence of the 
best-fit pseudo-continuum and Fe lines. The final measure 
of the pseudo-continuum flux, f\ (3000), provides the contin- 
uum luminosity L3000 which is further used to calculate Mbh 
and L/LEdd- It is important to note that the L3000 thus obtained 
is not the intrinsic, underlying AGN continuum since it also 
includes a contribution from the Balmer continuum emission, 
that is not accounted for by the fitting procedure. 

Once the pseudo-continuum and Fe contributions are sub- 
tracted, we fit the Mg II line itself. The model for the line 
consists of 3 Gaussians, two broad and one narrow compo- 
nent for each of the Mg II doublet lines. The intensity ra- 
tio of the doublet components is fixed to 1:1, suitable for 
optically-thick lines. The broad components are limited to 
line widths in the range 1,200 < FWHM < 10,000kms-', 
while the width of the narrow components is bound to 300 < 
FWHM < l,200kms"'. Although the relative contribution 
of the narrow components is very small (< 10%), the pro- 
cedure reproduces well the NLR width. This was verified 
by comparing the NLR FWHM resulting from the (separate) 
Mg II and H/3 fitting for the test sample of ^ 5,000 sources 



at 0.5 , 



, 0.75 mentioned above. We have also included a 



single absorption feature (also a Gaussian) to account for the 
frequent appearance of blue-shifted Mg II troughs. The best- 
fit model is used to measure the FWHM of the total broad 
single Mg II line and the total luminosity of both Mg II com- 
ponents. We stress that in earlier studies, the line width was 
calculated using the combined doublet profile, which has im- 
plications to the deduced Mbh (see § ^5. lb . The best-fit con- 
tinuum and Mg II line parameters are given in Table |2] The 
uncertainties that also appear in Table |2jreflect the true uncer- 
tainties associated with the spectroscopic reduction and line- 
fitting procedure. These were estimated from varying several 
of the parameters involved in these processes^ and examining 
the range of possible outcomes. We consider these uncertain- 
ties to be much more realistic than those derived solely from 
the flux noise. 

3.2. Estimating Mqh and L/L^dd 
The Lboi and Mbh estimates used in this work are based 
on the measured L3000 and FWHM(Mg ll). For the latter we 
adopt the McLure & Dunlop (2004) relationship: 



Mbh = 3.2x10^ 



^^3000 



T 0.62 r 



1044 



ergss 



FWHM (Mgll) 
lO^kms"' 



Mr 



(1) 

To estimate Lbob we have to assign a bolometric correction 
factor, /boi (A) = Lboi/ ALa- For this we adopt the Marconi et 
al. (2004) luminosity-dependent SED, which also provides a 
polynomial prescription for estimating /boi (4400). The pre- 

^ For example, varying flux calibrations for sources which were observed 
during several nights, the shape and normalization of the Fe II template and 
the possible identification of absorption features. 



scription can be used with other UV-optical monochromatic 
luminosities, by adopting the Vanden Berk et al. (2001) UV 
continuum (/^ i^"*'-^^). Our procedure relies on our own 
empirically-calibrated /boi (3000) vs. L3000 relation, derived 
using the 0.5 < z < 0.75 type-I SDSS AGN test sample. For 
each source, we measured both L5100 and L3000 and converted 
L5100 to Lboi using the relation of Marconi et al. (2004) and the 
Vanden Berk et al. (2001) template. This provides, for each of 
the ^ 5,000 sources, both L3000 and /boi (3000) — Lboi/^3000- 
The best-fit relation is, 

/boi (3000) = -0.58^3^000 44+ 3-85j2fooo 44- 8.38^3000.44+ 9.34 , 

(2) 

where ^3000,44 = log (L3000/ 10^^ ergss The bolometric 
corrections for our z — 4.8 sample range between 3.35 for 
the least luminous source and 3.43 for the most luminous 
one. The typical /boi (3000) is lower than those used in other 
studies (e.g., Elvis et al. 1994; Richards et al. 2006b) by a 
factor of about 1.5. This is mostly due to the fact that the 
Marconi et al. (2004) estimates of Lboi do not include most 
of the mid- to far-IR emission, which originates from dust 
around the central source. The normalized accretion rate is 
L/LEdd — Lboi/ (1-5 X 1O^**Mbh/M0) that is appropriate for 
ionized gas with solar metallicity. The deduced Mbh and 
L/LEdd are given in Table |2] 

We have also estimated Mbh from the C IV A 1549 line, that 
is observed in the SDSS spectrum of each source. The line 
fitting procedure was similar to the one described above for 
Mg II and also to those discussed in other studies (e.g., Shen 
et al. 2008; Fine et al. 2010). We calculated Mbh from L1450 
and FWHM(C IV) following the prescription of Vestergaard 
& Peterson (2006). We find that the C IV fine is systemat- 
ically broader than Mg II, similarly to the recent finding of 
Fine et al. (2010). The differences in widths translate to higher 
C iv-based Mbh estimates, with respect to those derived from 
Mg II. These findings confirm the suspicion about the prob- 
lematic use of the C iv-based method (see |T]i. A full analysis 
of these issues is deferred to a future publication. In what fol- 
lows, we rely solely on the Mg ll-based measurements of Mbh 
and L/LEdd- 

4. RESULTS 

4.1. Mbh and L/LEdd distributions 

The distributions of Mbh and L/LEdd for the z — 4.8 sam- 
ple are shown in Figure |2] They cover the range 8.03 < 
log(MBH/M0) < 9.82 and 0.18 < L/Lndd < 3.92. The me- 
dian values and 68%-percentiles (taken to be symmetric 
around the medians) correspond to (log (Mbh/M©)} = 8.89 

-034 and (L/LEdd)_= 0.59 JO ■ 

As explained in g3] both Mbh and L/LEdd are derived di- 
rectly from L3000, and both distributions are affected by the 
flux limit of the sample. The 1450A flux limit corresponds 
to /a (3000) > 2.1 X 10-'**ergscm-2s-i-i which, at z ~ 4.8, 
gives a limiting luminosity of L3000 — lO^^ergss"^. The lim- 
iting Mbh is estimated by combining the limiting L3000 and 
FWHM = 1, 500km s"' using Eq. [U This line width is the 
minimal width of any emission line for which the SDSS 
pipeline labels the target as a spectroscopic target of type 
"QSO"; a pre -requirement of our sample selection. We thus 
obtain log (Mbh/M0)jjj^J( ~ 8. Indeed, one of our sources 
(J1447-I-1025) has very similar properties, in terms of L300O1 
FWHM(Mg II) and Mbh- However, this source looks as an 
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Table 2 

Observed and Derived Properties 



Object ID (SDSS J) 


ZSDSS 


ZMgir' 


logi 


'1450*' 


log L3000 


L-qual. 


log ibol 


FWHM(Mg II) 


FWHM-qual. 


log Mbh 


log L/Lem 








(ergs 




(ergss"' ) 




(ergs s"' ) 


(km s"' ) 








J000749. 17+004119.4 


4.837 


4.786 


46, 


.39 


46.04 


3 


46.57 


3665 


3 


8.90 


-0.51 


J003525.28+004002.8 


4.757 


4.759 


46, 


.35 


46.38 


3 


46.91 


1805 


3 


8.49 


0.24 


J021043. 15-001818.2 


4.733 


4.713 


46, 


.66 


46.04 


2 


46.56 


4583 


3 


9.09 


-0.70 


J0331 19.67-074143.1 


4.738 


4.729 


46, 


.76 


46.55 


1 


47.09 


2360 


1 


8.83 


0.08 


J075907.58+1 80054.7 


4.861 


4.804 


46, 


.46 


46.54 


1 


47.07 


2717 


1 


8.95 


-0.05 


J080023.03+305 100.0 


4.687 


4.677 


46, 


.82 


46.73 


1 


47.26 


1404 


1 


8.49 


0.59 


J080715. 12+132804.8 


4.874 


4.885 


46, 


.71 


46.53 


2 


47.07 


3837 


3 


9.24 


-0.35 


J083920.53+352457.6 


4.777 


4.795 


46, 


.70 


46.24 


2 


46.77 


1971 


2 


8.49 


0.11 


J085707.94+321032.0 


4.776 


4.801 


46, 


.94 


46.72 


2 


47.25 


2851 


2 


9.10 


-0.03 


J092303.53+024739.5 


4.660 


4.659 


46, 


.33 


46.14 


1 


46.67 


2636 


1 


8.68 


-0.18 


J093508.50+080114.5 


4.699 


4.671 


46, 


.62 


46.33 


2 


46.87 


2714 


2 


8.82 


-0.13 


J093523.32+41 1518.7 


4.836 


4.802 


46, 


.66 


46.58 


2 


47.12 


3447 


3 


9.18 


-0.24 


J094409.52+100656.7 


4.748 


4.771 


46, 


.64 


46.40 


2 


46.93 


2128 


2 


8.65 


0.11 


J101759.64+032740.0 


4.917 


4.943 


46, 


.27 


46.10 


1 


46.63 


2822 


1 


8.71 


-0.26 


J105919.22+023428.8 


4.735 


4.789 


46, 


.65 


46.36 


2 


46.89 


2899 


2 


8.89 


-0.18 


Jll 1358.32+025333.6 


4.882 


4.870 


46, 


.49 


46.35 


2 


46.89 


3793 


3 


9.12 


-0.41 


Jl 14448.54+055709.8 


4.793 


4.790 


46, 


.13 


46.11 


3 


46.63 


3215 


2 


8.83 


-0.37 


J115 158.25+030341. 7 


4.698 


4.687 


46, 


.05 


45.91 


3 


46.44 


3741 


3 


8.84 


-0.57 


J120256.44+072038.9 


4.785 


4.810 


46, 


.28 


46.27 


1 


46.80 


2171 


1 


8.59 


0.04 


J123503.04-000331.6 


4.723 


4.700 


46, 


.07 


46.12 


2 


46.65 


4422 


3 


9.11 


-0.64 


J130619.38+023658.9 


4.852 


4.860 


46, 


.57 


46.82 


1 


47.35 


5340 


1 


9.71 


-0.54 


J131737.28+110533.1 


4.810 


4.744 


46, 


.39 


46.34 


2 


46.87 


3144 


2 


8.95 


-0.25 


J132110.82+003821.7 


4.716 


4.726 


46, 


.47 


46.17 


2 


46.70 


3651 


3 


8.98 


-0.45 


J132853.67-022441.7 


4.695 


4.658 


46, 


.42 


46.28 


2 


46.81 


3815 


2 


9.08 


-0.45 


J133125.57+025535.6 


4.737 


4.762 


46, 


.15 


46.02 


3 


46.55 


3445 


3 


8.83 


-0.46 


J134134.20+014157.8 


4.670 


4.689 


46, 


.87 


46.73 


2 


47.26 


6480 


2 


9.82 


-0.74 


J134546.97-015940.3 


4.714 


4.728 


46, 


.62 


46.08 


2 


46.60 


3412 


2 


8.86 


-0.43 


J140404.64+03 1404.0 


4.870 


4.903 


46, 


.55 


46.49 


2 


47.02 


5360 


2 


9.51 


-0.66 


J143352.2 1+0227 14.1 


4.721 


4.722 


47, 


.14 


46.84 


2 


47.37 


2622 


2 


9.11 


0.09 


JI43629.94+063508.0 


4.850 


4.817 


46, 


.62 


46.44 


2 


46.98 


3052 


2 


8.99 


-0.19 


J144352.95+060533.1 


4.879 


4.884 


46, 


.44 


46.16 


3 


46.69 


3609 


3 


8.96 


-0.45 


J144734. 10+1025 13.2 


4.686 


4.679 


46, 


.29 


45.99 


1 


46.51 


1407 


1 


8.03 


0.30 


J151 155.98+040803.0 


4.686 


4.670 


46, 


.62 


46.32 


2 


46.86 


1735 


2 


8.42 


0.26 


J161622.1 1+050127.7 


4.872 


4.869 


47, 


.08 


46.80 


2 


47.33 


3910 


1 


9.43 


-0.27 


J165436.86+222733.7 


4.678 


4.717 


47, 


.14 


46.48 


1 


47.02 


5637 


1 


9.55 


-0.70 


J205724.15-003018.0 


4.663 


4.680 


47, 


.04 


46.83 


2 


47.36 


3050 


2 


9.23 


-0.05 


J220008.66+001744.8 


4.818 


4.804 


46, 


.70 


46.51 


2 


47.04 


2404 


2 


8.82 


0.04 


J221705.72-001307.7 


4.689 


4.676 


46, 


.44 


46.28 


3 


46.81 


2274 


3 


8.63 


0.00 


J222509. 16-001406.8 


4.888 


4.890 


46, 


.97 


46.70 


1 


47.23 


3504 


1 


9.27 


-0.21 


J224453.06+134631.8 


4.657 


4.656 


46, 


.30 


46.06 


2 


46.58 


2516 


2 


8.58 


-0.17 



" Redshift measured from the best-fit models of the Mg II lines. 
Monochromatic luminosity at rest-wavelength 1450A, obtained from the SDSS/DR6 spectra and redshifts. 

Quality flag associated with Ljooo- Quality flags of "1", "2" and "3" correspond to calibration and/or continuum measurement uncertainties of ~ 10%, '■ 
and ~ 40%, respectively. 

Quahty flag associated with FWHM(Mg II). Quality flags of "1", "2" and "3" correspond to uncertainties of ~ 10%, 20—30% and ~ 50%, respectively. 



20% 



outlier in the Mbh distribution. The next smallest source has 
log (Mbh /M0) — 8.4, which represents the low end of the 
Mbh distribution much better. We conclude that our sample is 
complete down to log (Mbh/M©) ~ 8, and the fact that we ob- 
serve (almost) no sources with log (Mbh/M©) < 8.5 is a real 
characteristic of the population of spectroscopically observed, 
optically selected type-I AGNs at z — 4.8. 

The limiting L/Ledd is deduced by combining the lowest 
L3000 with FWHM (Mgll) = 10,000kms"\ the upper bound- 
ary tha t FW HM(Mg 11) can take in our line-fitting procedure 
(see § ^3.11 ). The resulting limiting value is (L/LEdd)iimit — 
0.04. Clearly, the limiting L/Ledd is a factor of ^ 4.5 below 
the lowest L/Lndd we actually observe. This reflects the fact 
that the broadest Mg II Une we observe has FWHM (Mgll) ~ 
6,500kms-^ 

The highest observable L/LEdd is determined by a com- 
bination of the narrowest observable FWHM(Mg ll) and the 
highest luminosities. For instance, the faint and narrow- 
lined source J1447H-1025 is one of the higher L/Lem sources 
in our sample. On the other hand, the hightest L/L^Ad 



source (J0800H-3051) has a similarly narrow Mgll line but 
is among the most luminous sources in the sample (a factor 
of ~ 5.6 more luminous than J1447H-1025), thus approaching 
L/Leaa — 4. Since there are no brighter z — 4.8 sources in the 
entire SDSS/DR6 database, we conclude that our sample un- 
veils the highest-L/Ledd (—4) sources atz ~4.8. We also note 
that ^ 1/4 of the AGNs in our sample have 1 < L/Lndd ^ 4. 
Several studies (e.g., Mineshige et al. 2000; Wang & Netzer 
2003; Kurosawa & Proga 2009) have shown that such high 
L/Leaa can be produced in several types of accretion disks. 



4.2. The growth times of SMBHs from z = 20 to z — 4.8 

Given the observed Mbh, i'boi ™d L/Leaa, we can con- 
strain the lifetimes of the SMBHs in our sample. We first 
assume that accretion proceeds with a constant L/Leaa- This 
results in an exponential growth of Mbh, with an e-folding 
time ("Salpeter time") of 



T = 4x 10' 



yr, 



(3) 
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log^Bo/LEdd 

Figure 2. Distributions of Mbh (top) and L/Ledd (bottom) for the z — 4.8 
sample. 

where 77 is the radiative efficiency (Lboi = ?7Mnfaiic^). The 
time required to grow by such accretion from a seed BH of 
mass Mseed to Mbh is 



: T In yr. 



(4) 



Despite several attempts to estimate ry (e.g., Volonteri et 
al. 2005; King, Pringle & Hofmann 2008, and references 
therein), its typical value and redshift dependence are poorly 
constrained. A broad range of 77 ~ 0.05-0.3 is required to 
match the local BH mass function to the integrated accretion 
history of SMBHs (e.g., Shankai-et al. 2004; Wang et al. 2009; 
Shankar et al. 2010a). There is an even broader range of pos- 
sible values of Mseed, depending on the various mechanisms to 
produce such objects. Remnants of population-Ill stars would 
result in Mseed ^ 10- lOOMg (e.g. Heger & Woosley 2002), 
while collapse models of either dense stellar clusters or gas 
halos predict seed masses as large as Mseed ^ 10^- lO^'M© 
(e.g., Begelman, Volonteri & Rees, 2006; Devecchi & Volon- 
teri 2009; see Volonteri 2010 for a review). For the sake of 
consistency with earlier z > 2 studies (N07), we assume here 
r/ = 0.1 and Mseed — lO'^M©. With this choice of 77, the e- 
folding times in our sample range from ^ 10 to 240 Myr 
and the faster-accreting 50% of the sources show t ^ 75 Myr. 
Such SMBHs can increase their mass by as much as a factor 
of ~ 790 within ^ 500 Myr. We compared the growth of the 
SMBHs in our sample up to z ~ 4.8 with the corresponding 
age of the Universe at that redshift (funiverse — 1-21 Gyr). The 
assumption of continuous, constant L/Ledd growth for such 
a long period is, of course, somewhat simplistic, given the 
timescales of standard accretion disks (a review of this and 
related issues can be found in King 2008). About 65% of 
the sources have f growth < ^Universe- This is a larger fraction 
than observed in lower-redshift samples, which suggests that 
z ~ 4.8 might represent the sought after episode of fast growth 
for most SMBHs. 

We also calculated the regress of Mbh with increasing red- 
shift assuming the aforementioned exponential growth sce- 




t„ ■ [Gyr] 

Universe ^ 

Figures. The calculated growth of Mbh from z = 20 to z — 4.8. Black 
squares represent the masses of the z — 4.8 sources and the dotted lines show 
the regress in Mbh, assuming exponential growth with the observed L/Lem 
and r; = 0. 1 . The thick lines illustrate a choice of different values of rj for a 
particular source (J1436+0635) which has Mbh and L/Leu that resemble the 
sample's median values. 

nario, up to z = 20. These evolutionary tracks are shown 
in Figure [3] The point at which such tracks cross the y- 
axis in Fig [5] can be regarded as the value of Mseed required 
to match the observed Mbh at z — 4.8, given the assump- 
tion of a constant L/Ledd growth. The diagram demonstrates 
how the z — 4.8 sample can constrain the different BH for- 
mation and growth scenarios. For example, about 40% of the 
z ~ 4.8 sources could have grown from seed BHs with Mseed < 
IOOM0, i.e. stellar-remnant BHs. A further ^ 20% could 
have grown from BH seeds with 1,000 < Mseed < 10^ M©, and 
another ^ 20% from BH seeds with 10^ < Mseed < lO^M©. 
The remaining, low est- L/L^^a sources (^ 10%) have 10^ < 
Mseed < 10** Mf;,|3. Since even the most extreme scenarios can- 
not produce such high Mseed (Volonteri 2010), these sources 
had to accrete at higher rates in the past in order to attain 
their measured mass. Another explanation is a combination 
of massive seed BHs with low spins (i.e. small rj) at very 
early epochs. For example, if we assume rj = 0.05 then at 
z = 20 all the sources have an implied Mseed < lO^M©. As 
mentioned above, the radiative efficiency could be higher, e.g. 
^ 0.2-0.3. Assuming such high values of 77 in our calcu- 
lations naturally leads to the requirement of more massive 
seed BHs. For example, in the extreme case of 77 = 0.3 in 
all sources, we find that 95% of the seed BHs would have 

Mseed > 5OOM0. 

5. DISCUSSION 

5.1. Comparison with other studies 

Our main goal is to identify the epoch at which most 
SMBHs experienced the first episode of fast growth. We focus 
on a comparison of our z — 4.8 sample to several z > 2 sam- 
ples with reliably measured Mbh and L/LEdj. In S04 and N07 

^ Alternatively, assuming Mj^ed = lO'^M.,, these SMBHs have growth 
times which are about twice the age of the Universe at the corresponding 
redshift. 
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Figure 4. Ltoi, Mbh and L/L^dd vs. redshift, for samples of different red- 
shifts discussed in the text: the new z. ~ 4.8 sample presented here (black 
squares), the z — 2.4 and ~ 3.3 samples of S04 and N07 (magenta and blue 
triangles) and the combined z — 6.2 sample from K07 and WIO (red circles). 

we studied a sample of 44 z ~ 2.4 and ~ 3.3 type-I AGNs us- 
ing NIR spectroscopy to measure L5100 and FWHM(H/3). The 
study showed that a large fraction of the high Mbh sources ac- 
crete at a rate that is well below the Eddington limit. Assum- 
ing these sources accrete at constant L/L^^i, such accretion 
rates cannot explain their measured masses and are in contra- 
diction with several theoretical predictions. Having obtained 
new data on the z — 4.8 sample, we can now compare several 
groups of AGNs, at several epochs, in an attempt to follow 
their growth all the way from z ~ 6.2 to z ~ 2.4. 

The samples we consider here are the new one at z ~ 4.8 
and our earlier (S04 & N07) samples at z — 2.4 and ~ 3.3. 
We also consider a small number of sources at z — 6.2 from 
the samples studied by K07 and WIO, which have 5 and 9 
reliable Mbh & L/Leaa measurements, respectively. These 
small samples do not have a common flux limit and simply 
represent the up-to-date collection of the sources discovered 
at those redshifts that were also observed in follow-up NIR 
spectroscopy. As such, they are not representative of the AGN 
population at z — 6.2. This, along with the small size of the 
samples, limits their statistical usefulness. To use the data 
from the z — 6.2 samples, we re-calculated Mbh and i'boi us- 
ing the methods described above. First, we corrected the pub- 
lished FWHM(Mg II) according to the assumption of only one 
of the doublet components. This reduces FWHM(Mg ll) by a 
mean factor of ^ 1 .2 and the reported Mbh by a mean fac- 
tor of 1 .44. Second, the smaller /boi (3000) adopted here 
reduces the reported Lboi by a factor of ^ 1.49. Regarding 
L/Lgdd^ the two corrections almost completely cancel out. 

In Figure |4] we present Lbob Mbh and L/LEdd for the new 
z ~ 4.8 sample, the z — 2.4 and ~ 3.3 samples of S04 and 
N07, and the combined sample of z — 6.2 AGNs. The four 
samples suggest an increasing Mbh and decreasingL/LEdd 
with cosmic time. This is better illustrated in Figure^ where 
we present the cumulative distribution functions (CDFs) of 
Mbh and L/Ledd- The increase in the median Mbh and the 
decrease in the median L/L^dd as a function of redshift are 




Figure 5. The cumulative distribution functions (CDFs) of Mbh (top) and 
L/LEdd (bottom) for the samples discussed in the text. 

evident. In particular, there is a clear shift of ^ 0.5 dex be- 
tween the median Mbh values of the z — 4.8 and z — 3.3 sam- 
ples, in the sense of a lower Mbh at z — 4.8. There is also an 
opposite shift between the median L/LEdd values. The differ- 
ences between the z — 3.3 and the z — 2.4 samples are much 
smaller, although the z — 2.4 sample includes several sources 
with extremely massive BHs (log (Mbh/Mq) > 10) which are 
not observed at earlier epochs. Regarding Mbh, only ^ 14% 
of the z — 3.3 sample (2 sources) lie below the median value 
of the z ~ 4.8 sample (log (Mbh/Mq) = 8.92). Similarly, only 
^ 13% of the z — 4.8 sample (5 sources) lie above the median 
value of the z — 3.3 sample (log (Mbh/Mp^) = 9.37) and even 
less (^ 6%; 2 sources) above the median value of the z — 2.4 
sample (log(MBH/M0) = 9.59). 

We have further tested the significance of these differences 
by performing a series of two-sample Kolmogorov-Smirnov 
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tests. The null hypothesis that the observed distributions of 
L/Leaa at z ~ 4.8 and at z ~ 3.3 (or at z ~ 2.4) are drawn 
from the same parent distribution is rejected with significance 
levels > 99%. A similar test was applied to the z — 4.8 and 
z ~ 6.2 samples, and could not reject the null hypothesis, sug- 
gesting that the distributions of L/Ledd at z — 4.8 and z — 6.2 
are statistically similar. We obtain similar results when com- 
paring the distribution of Mbh at z — 4.8 to that of the z — 2.4 
and 2± 3.3 samples. Due to the large fractionof z — 6.2 sources 
with log (Mbh/-Wp) ^ 8.5, we can reject the hypothesis that 
the Mbh values at z — 4.8 and z — 6.2 are drawn from the same 
distribution. 

Some of the above results may be biased by the fact that the 
four samples cover a different range of Lboh which originate 
from the different target selection criteria. We thus repeated 
the statistical tests aforementioned, focusing on subsamples 
which share a common range of 46.4 < log (Lboi/ergss"') < 
47.4 (i.e., matching the Lboi range of the z — 4.8 sample). All 
but one comparison result in the same conclusions, with sim- 
ilar confidence levels. The one exception are the distributions 
of Mbh at z — 4.8 and z — 6.2, for which the null hypothe- 
sis now cannot be rejected, suggesting that these distributions 
represent the same parent population. This is not surprising 
given the lower luminosities of most of the z — 6.2 sources, 
the dependence of Mbh on source luminosity, and the incom- 
pleteness of the z — 6.2 sample. 

We conclude that there is strong evidence for a rise in Mbh 
and a drop in L/LEdd of about a factor of 2.8 between z — 4.8 
and z — 3.3. This strong trend is not observed with respect to 
neither lower nor higher-redshift samples, which span similar 
periods of time 0. It thus seems that the most massive BHs, 
associated with the most luminous AGNs, started the episode 
of fast BH growth at redshifts above about z ^ 5. By z ~ 
2 — 3, these SMBHs reach their peak (final) mass ( IO^OMq) 
and their mass accretion is less efficient. 

5.2. The growth of active SMBHs from z — 4.8 to z — 3.3 and 

Given the trends and differences in Mbh and L /L^ dd, as well 
as the short e-folding times of the z ~ 4.8 AGNs (§ ^4.21 1. it is 
possible to think of the z — 4.8, z — 3.3 and z — 2.4 samples as 
representing different evolutionary stages of the same parent 
population of SMBHs. In what follows, we focus on the z — 
4.8 and z — 3.3 samples (separated by ^ 680Myr) to test this 
evolutionary interpretation, assuming various scenarios. 

We consider two growth scenarios: constant accretion rate 
(i.e. constant Lboi) and constant L/LEdJZl As explained, 
the assumption of constant L/L^dd results in an exponential 
growth of Mbh and given our chosen value of 77, the mass 
growth e-folding time is r ~ 45 {L/L^ddY^ Myr, which trans- 
lates to T < 240 Myr (L/Lndd > 0.18) for the z ~ 4.8 sources. 
Thus, even the lowest L/Ledd SMBHs at z — 4.8 could have 
increased their Mbh by a factor of 20 by z — 3.3, which is 
much larger than the typical difference in Mbh between the 
two samples. The assumption of a constant Lboi results in 
much slower growth. The luminosities of the z — 4.8 sources 
translate to 4 ^Mbh ^ 37M0yr~^ Even this much slower 
growth scenario results in very high Mbh at z — 3.3. This sce- 

* For the adopted cosmology, the physical time between z = 4.8 and z = 3.3 
is ~ 680Myr and between z = 3.3 and z = 2.4 is ^ 790 Myr. 

' There are, of course, many other possible scenarios involving, for exam- 
ple, host related evolution. These are beyond the scope of the present paper. 
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Figure 6. Evolution scenaiios for the z — 4.8 SMBHs. Symbols are identical 
to those in Fig. |4] Solid lines describe the growth of Mbh under the as- 
sumption of a constant Lboi, while dotted lines represent the constant L/Ledd 
scenario. 

nario also produces very low accretion rates (L/Ledd < 0. 1) at 
z~3.3. 

In Figure |6] we present evolutionary tracks for our z — 4.8 
AGNs to z — 3.3 and eventually to z = 2. Clearly, continu- 
ous constant Lboi growth results in too large Mbh and can- 
not reproduce the lower Mbh sources at z — 2.4 and ~ 3.3. 
Specifically, the calculated distribution of Mbh for the z — 4.8 
sources, when evolved to z — 3.3, has a median value of 
log (Mbh/Mtj) ~ 10, larger by ^ 0.6 dex. than the observed 
median of the z — 3.3 sample. This scenario also fails to repro- 
duce the observed range of L/Ledd, since the lowest L/Ledd 
sources at z — 3.3 have L/Ledd ^0.1. As the dotted evolu- 
tionary tracks in Fig |6] demonstrate, the constant L/Lndd sce- 
nario produces even larger masses, and is only feasible if all 
the z — 4.8 sources cease their accretion shortly after their ob- 
served active phase. This scenario assumes a constant L/Lgdd^ 
hence the significant difference between the L/Ledd distribu- 
tions at z — 4.8 and at z — 3.3 (^ 0.5 dex ; see Fig|5]) is not 
resolved. We also note that in both scenarios the z — 4.8 
sources would have been easily observed at z — 3.3, since 
they would have Lboi which is either similar (linear growth) or 
much higher than (exponential growth) the sources observed 
by S04 and NO? at z — 2.4 and ~ 3.3. From these two sim- 
plistic growth scenarios we conclude that the fast growth of 
the z — 4.8 sources must either experience a shut down before 
z ~ 3.3, or accrete in several short episodes, with duty cycles 
which are much smaller than unity. 

To further constrain the evolution of the observed SMBHs, 
we ran a series of calculations with different duty cycles for 
each of the two evolutionary scenarios. In each calculation 
we assembled the distribution of calculated Mbh at exactly 
z = 3.3 and z = 2.4. Several of these distributions are shown 
in Figure |7] For the fast, constant L/Lgdd scenario, the only 
calculated distributions which resemble the observed distri- 
bution of Mbh at z — 3.3 are those with a duty cycles in the 
range 7.5- 12.5%. The observed distribution at z — 2.4, on the 
other hand, can only be achieved by duty cycles of 5-7.5%. 
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Figure 7. Observed and calculated Mbh CDFs for different evolutionary sce- 
narios. In all panels solid black lines show the observed Mbh CDF at z ~ 4.8, 
while magenta lines (left panels) and blue lines (right panels) show the CDFs 
at z ~ 3.3 and at z — 2.4, respectively. Dashed black lines are the calculated 
CDF of the z — 4.8 sample assuming different growth scenarios and duty cy- 
cles. Top panels assume a constant L(,oi scenario, and duty cycles of 5, 20, 
and 50% (from left to right). Bottom panels assume a constant L/Z-Edd sce- 
nario, and duty cycles of 5, 10 and 15% (from left to right). The dotted lines 
in the bottom panels illustrate how a duty cycle of 50% produces extremely 
over-massive BHs at z — 3.3 and z — 2.4. All the calculations assume 77 = 0.1. 
Assuming the extreme cases of = 0.05 (or 0.3) would mean that the plotted 
calculated CDFs correspond to duty cycles which are a factor of ~ 2 smaller 
(or larger; see text for details). 



For the constant Lboi scenario, the z ~ 3.3 distribution can be 
matched by assuming a duty cycle in the range of 15-25%, 
while the z — 2.4 distributions can be partially explained by 
assuming duty cycles of ^5-25%. The ranges of duty cy- 
cles which account for the observed distribution of L/Ledd are 
somewhat different: 10-15% for the L/L^^i distribution at 
z ~ 3.3 and 5 - 10% for the one at z ~ 2.4. All the above cal- 
culations assumed rj — Q.l. Naturally, lower (higher) radiative 
efficiencies will require shorter (longer) duty cycles, to repro- 
duce the same calculated CDFs at z ~ 2.4 and ~ 3.3. This 
means that for any assumed 77 in the range 0.05 < ry < 0.3 the 
aforementioned "acceptable" duty cycles can change by up to 
a factor of ^ 2, where the exact factor scales as ?;/ (1 -?7). We 
note that Fig. |7] suggests that, in both evolution scenarios, the 
more massive BHs at z — 2.4 and ~ 3.3 seem to have grown at 
higher duty cycles than the less massive ones. Alternatively, 
this can be interpreted as an increase in radiative efficiency 
with increasing resultant Mbh at z — 2.4 and ~ 3.3. In partic- 
ular, the observed distributions of Mbh and L/Lem at z — 2.4 
are much more complex than the calculated ones and no sin- 
gle, fixed duty cycle can account for the shape of the observed 
distributions. The reason for this apparent discrepancy might 
also be the way the z — 2.4 sources were selected in the S04 
and N07 studies. 

To illustrate how the above duty cycles facilitate an evolu- 
tionary connection between the three samples, we present in 
Figure [8] the evolutionary tracks of the z — 4.8 sample, sim- 
ilar to Fig|6] but this time assuming duty cycles of 10% and 
20%, for the constant L/Ledd and constant Lboi scenarios, re- 
spectively and assuming again 77 = 0.1. We also verified that 
evolving the Mbh of the z — 2.4 and ~ 3.3 sources backwards. 
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Figure 8. Evolution scenarios for the z — 4.8 SMBHs with various duty cy- 
cles. Symbols are identical to those in Fig. |4] Solid lines describe Mbh 
growth under the assumption of constant Lboi and a duty cycle of 20%. Dot- 
ted lines represent the constant L/Ledd scenario and a duty cycle of 10%. 

under the assumption of constant L/Lem and a duty cycle of 
10% results in a distribution of Mbh which is not critically 
different than the one directly observed at z — 4.8. However, 
we again find that the few highest Mbh sources at z — 2.4 
probably require higher duty cycles. 

We conclude that all the observed measurements of Mbh, 
L/Lem and Lboi are consistent with duty cycles of about 10- 
20%. Considerably longer duty cycles can be consistent with 
observations only if we assume an extremely high radiative 
efficiency (rj ~ 0.3) for all the z — 4.8 sources. Such low duty 
cycles are in good agreement with the models presented in 
Shankar et al. (2009), which are able to reproduce the growth 
of a similar population of SMBHs, in terms of redshift and 
range of Mbh (cf. their Fig. 7).. Our constrains on duty cy- 
cles are, however, in contrast with several other studies that 
suggest that the duty cycle at high redshift should reach unity, 
based on the clustering of high-redshift AGNs (e.g. White et 
al. 2008; Wyithe & Loeb 2009; BonoH et al. 2010, and refer- 
ences therein). 

All this leads us to suggest that for the most massive 
z > 2 type-I AGNs, an epoch of fast SMBH growth took 
place before z ~ 4, which we partially observe at z — 4.8. 
This epoch is efficient enough to produce the very massive 
(log(MBH/M0) > 10) BHs observed at z ~ 2.4. The fast 
growth slows down before z — 3.3, and the sources observed 
at z — 3.3 show much lower L/Ledd- In addition, there is no 
significant rise in Mbh between z — 3.3 and z — 2.4. The 
growth of Mbh during this epoch can only be explained by as- 
suming that mass accretion proceeded in short episodes, last- 
ing an order 10-20% of the total period. This translates to 
accretion episodes which (cumulatively) last ^ 70- 140 Myr ( 
^ 150- 300 Myr) over the period between z — 4.8 andz — 3.3 
(z — 4.8 and z — 2.4), respectively. Another possibility is 
that the z — 4.8 SMBHs have experienced an even faster shut 
down, and so at z — 3.3 their inactive relics have masses which 
do not differ significantly from those we observe at z — 4.8. 

Major mergers between massive, gas-rich galaxies are ca- 
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pable of supplying large amounts of cold gas directly to the 
innermost regions, to be accreted by the central SMBHs. De- 
tailed simulations of such mergers suggest that these events 
can fuel significant BH growth over periods of ^ 1 Gyr. The 
mass accretion rate and source luminosity during such merg- 
ers may vary on time-scales of few x lOMyr (Di Matteo et al. 
2005; Hopkins et al. 2006; Hopkins & Hernquist 2009). Ac- 
cording to Hopkins et al. (2006), AGNs would appear to have 
Lboi ~ 2.7 X 10^^ ergs s~' (i.e., matching the range we observe 
at z ~ 4.8) for typically ^ 100 Myr. Numerical studies sug- 
gest that massive dark matter halos may undergo more than 
one major mergers per Gyr at z ~ 4.8 (e.g. Genel et al. 2009 
and references therein). Thus, during the ^ 680 Myr between 
z ~ 4.8 and z ~ 3.3, the hosts of the z ~ 4.8 SMBHs may have 
undergone a single merger, during which the period of signifi- 
cant accretion by the SMBHs would last for ^ 100 Myr. This 
is in good agreement with the duty cycles of 10-20% we find 
here, which correspond to 70- 140 Myr. The decline in activ- 
ity observed for the most massive BHs at z — 3.3 and z — 2.4 
may then be associated with the decline in the major merger 
rate, which drops by a factor of ^ 4 between z — 4.8 and 
z ~ 2.4 (Genel et al. 2009). However, it is likely that the hosts 
of the z — 2.4 SMBHs have experienced an additional major 
merger during the ^ 790 Myr between z — 3.3 and z — 2.4, 
especially if these SMBHs reside in the more massive dark 
matter halos. In such a scenario, the largest BHs at z — 2.4 
may have gathered their high mass during more efficient (i.e., 
higher duty-cycle) accretion episodes. This requirement of 
higher duty-cycles for higher Mbh sources is also reflected in 
our analysis (see upper-right panel of Fig. |7]). If major merg- 
ers are indeed the main drivers of SMBH accretion history at 
z ^ 2-5, the results presented here predict that the host galax- 
ies of the z — 4.8 sources would be found in the early stages of 
major mergers, while the hosts of z — 2.4 and ~ 3.3 sources 
would appear to be in either later stages, or in merging sys- 
tems which have a lower mass ratio. 

There are several other mechanisms to make large amounts 
of cold gas available for BH accretion. Most of these might 
also trigger intense star formation (SF) activity in the hosts 
of SMBHs. Indeed, many high-luminosity AGNs show ev- 
idence for intense SF (e.g. Netzer 2009b; Lutz et al. 2010, 
and references therein). New observations by Herschel sug- 
gest that SF in the hosts of the most luminous AGN peaks 
at z ^ 3 and quickly decreases at later epochs (Serjeant et al. 
2010). If correct, it will indicate that the amount of gas avail- 
able for both BH accretion and SF has depleted by z 2 - 3, 
consistent with the decrease in SMBH accretion activity we 
find here. Future Herschel and ALMA observations of the 
hosts of the z — 4.8 AGNs may be able to reveal the presence 
of such SF activity, and perhaps the amount and dynamical 
state of the cold gas. This will enable a better understanding 
of the galaxy-scale processes which drive the accretion his- 
tory of the fast-growing z ~ 4.8 SMBHs. 

6. SUMMARY 

We present new //-band spectroscopy for a flux-limited 
sample of 40 type-1 SDSS AGNs at z ~ 4.8. The sample 
covers ~ 1/4 of all the (spectroscopically observed) SDSS 
sources at that redshift band, and thus about ^ 1/20 of the to- 
tal population of Lboi > 2.75 x lO^^ergss"^ sources (over the 
entire sky). The main results of our study are: 

1. The z — 4.8 AGNs have, on average, higher accretion 
rates and lower massees than those observed at lower 



redshifts. The accretion rates and masses are compara- 
ble to those of the small, incomplete samples of z — 6.2 
AGNs. 

2. We have observed an epoch of fast SMBH growth, 
probably the very first such phase for most SMBHs. 
Assuming continuous growth from about z = 20, these 
observations provide the very first look at the distribu- 
tion of seed BHs in the early universe. About 65% of 
the SMBHs at z — 4.8 have had enough time to grow 
to their observed Mbh, assuming continuous accretion 
at the observed L/Lgcid- About 40% of the sources 
could have started their growth from BH seeds which 
are stellar remnants (Mjeed < lOOM©). For the minor- 
ity the sources, those with small L/LEdd, there might 
have been an ever earlier epoch of faster accretion or, 
perhaps, they started their growth from a much larger 
seeds (M^eed - 10"^ M0). 

3. The z — 4.8 sources can be regarded as the progenitor 
population of the most massive (Mbh ^ lO^^M©) BHs 
observed at z — 2.4 and ~ 3.3. The growth rate of those 
massive BHs seems to be much slower between z — 3.3 
and z — 2.4. 

4. A comparison of the observed distributions of Mbh at 
z ~ 4.8, z ~ 3.3 and z ~ 2.4 indicates that the z ~ 4.8 
sources either completely stop their accretion shortly 
after z — 4.8, or that their accretion proceeds in rel- 
atively short episodes. We find that for mass growth 
rates that follow either constant L/Lgdd or constant Lboi 
scenarios, duty cycles of either ~ 10% or ~ 20%, re- 
spectively, give reasonable agreement to the observed 
distiibutions of Mbh- 
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